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Abstract 



Projective space of order n over a finite field GF{q), denoted by 'Pq{n), is a set of all the subspaces 
of a vector space GF{q)'^. The projective space is a metric space with the distance function 
ds{U,V) = d\m{U) + dim(F) - 2dim(L'" fl V), for all U,V e Vq{n). A code in the projective 
space is a subset of Vq{n). Koetter and Kschischang showed that codes in projective space are 
useful for errors and erasures correction in random network coding. If the dimension of each 
codeword is restricted to a fixed integer, the code forms a subset of a finite-field Grassmannian. 
Such a code is called a constant-dimension code. Constant-dimension codes are analogous to 
constant-weight codes in the Hamming space. 

In our work we introduce a representation of codewords in Grassmannian by codewords in 
the associated constant-weight codes and matrices in reduced echelon form. We describe an 
algorithm for detecting the distance between two subspaces of Vq{n) based on such represen- 
tation. It is known that the size of Grassmannian can be determined by Gaussian coefficient. 
Using the connection between the Gaussian coefficient and the number of partitions we show 
encoding methods for constant-dimension codes, based on representation of partitions by Ferrer 
diagrams. 

Next, we consider rank-metric codes. A rank-metric code C is a subset of GF{(f^Y' with 
the rank-distance ci_R(ci,C2) = rank(ci — C2), for ci,C2 G GF{q'^)^. We construct new constant- 
dimension codes, based on optimal rank-metric codes. 

We continue by studying the codes in unrestricted projective space. We generalize the 
representation of codewords in unrestricted projective space by vectors in the Hamming space 
and matrices in reduced echelon form. We construct new codes in Vq{n), using our construction 
of codes in Grassmannian. 
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Chapter 1 

Introduction 



Let GF{q) be a finite field of size q and V a vector space of dimension n over GF{q). Pro- 
jective space (or projective geometry) of order n over GF{q), denoted by Vq{n), is a set of all 
subspaces of vector space V . The set of all /c- dimensional subspaces of F, for any /c < n is a 
Grassmannian, denoted by Qq{n,k). The size of Grassmannian is equal to the g-ary Gaussian 



coefficient, denoted by 



n 
k 



. A projective space Vq{V) can be viewed as a union of Grassman- 



q 

nians, i.e. Vq{n) = U Qq{n, k). The projective space is a metric space with distance function 
ds{U, W) = dim([/) + dim(VF) - 2dim(C/ n W), for all U,W ^ 'Pgin). Code in projective space 
C C Vq{n) of size M and minimum distance d is denoted by (n, M, d) code. If an (n, M, d) 
code C is contained in Qq{n,k) for some k, than C denoted by (n, M, d, k) code and called a 
constant-dimension code. Constant-dimension codes are analogous to constant-weight codes in 
the Hamming space. 

Our interest in codes in projective space is motivated by their applicability to errors and 
erasures correction in random network coding, as was shown by Koetter and Kschischang [1]. 
In their work they introduce an operator channel with input and output alphabet Vq{n) that 
can be viewed as a standard discrete memoryless channel, and prove that an (n, M, d) code C 
in projective space can correct any t errors and any p erasures if 2{t + p) < d. 

The construction of codes in projective space can be facilitated by exploiting their connection 
with rank-metric codes OlSjl^. A rank-metric code C is a subset of GF{q^)'^ with rank-distance 
dR{ci,C2) = rank(ci -C2), for ci,C2 G GF{q"^Y. 

The focus of this work is to analyze the codes in projective space. In the following we provide 
a high level overview of random network coding and rank-metric codes. 



1.1 Random network coding and operator channel 

A network is a directed graph, where the edges represent pathways for information. Using the 
max-flow min-cut theorem, one can calculate the maximum amount of information that can be 
pushed through this network between any two graph nodes. It was shown that simple forwarding 
of information between the nodes is not capable of attaining the max-flow value. Rather, by 
allowing mixing of data at intermediate network nodes this value can be achieved. Such encoding 
is referred to as network coding [5]. 
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In the example in Figure 11.11 two sources having access to bits A and B at a rate of one bit 
per unit time have to communicate these bits to two sinks so that both sinks receive both bits 
per unit time. All links have a capacity of one bit per unit time. The network problem can 
be satisfied with the transmissions specified in the example but cannot be satisfied with only 
forwarding of bits at intermediate packet nodes. 




1.1.1 Operator Channels and Coding for Operator Channels HI 

In random network coding, the communication between transmitter and receiver occurs in series 
of rounds or "generations"; during each generation, the transmitter injects a number of fixed- 
length packets into the network, each of which may be regarded as a vector of length over 
a finite field GF{q). These packets propagate through the network, possibly passing through 
a number of intermediate nodes between transmitter and receiver. Whenever an intermediate 
node has an opportunity to send a packet, it creates a random GF{q)- linear combination of 
packets it has available and transmits this random combination. Finally, the receiver collects 
such randomly generated packets and tries to infer the set of packets injected into the network. 
There is no assumption here that the network operates synchronously or without delay or that 
the network is acyclic. 

Let {pi,P2, •••)Pa/}, Pi £ GF{q)^ denote the set of injected vectors. In the error-free case, 
the receiver collects packets yj = Yli=i ^j-iPi^ ^ ^ j ^ L, with unknown, randomly chosen 
coefficients hj^i G GF{q). 

^This section is borrowed from [T]. 
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If we choose to consider the injection of T erroneous packets, this model is extended to 
include error packets ej, 1 <t <T to give 

M T 
i=l t=l 

where gj^t G GF{q) are unknown random coefHcients. 

The matrix form representation of the transmission model is 

y = Hp + Ge, 

where H and G are random L x M and L xT matrices, respectively, p is the M x N matrix 
whose rows are the transmitted vectors, and e is the T x N matrix whose rows are the error 
vectors. 

Let W he a fixed iV- dimensional vector space over GF{q), and let Vg{W) be the set of all 
subspaces of W, i.e., the projective space of order N. 

Definition 1.1.1. An operator channel C , associated with the ambient space is a channel 
with input and output alphabet VqiW). The channel input V and channel output U are related 
as 

U = Hk{V)®E, 

where T-Lk is an erasure operator, E G VqiW) is an arbitrary error space, and © denotes the 
direct sum. If dim(y) > k, the erasure operator Tik acts to project V onto a randomly chosen 
/c- dimensional subspace of V\ otherwise, Hk leaves V unchanged. If the erasure operator Hfe 
satisfies dim(y) — (\\m{l-ik{V)) = p, we say that Hjt corresponds to p erasures. The dimension 
of E is called the error norm t{E) of E. 

A code C for an operator channel with ambient space W is simply a nonempty subset of 
VqiyV). The minimum distance of C is denoted by 

d{C):= min dsiX.Y). 

X,YeC,Xj^Y 

The maximum dimension of the elements of C is denoted by 

1{C) := maxdim(X). 

The importance of the minimum distance d{C) for a code C C V(W) is given in the fol- 
lowing theorem, which provides the combined error-and-erasure-correction capability of C under 
minimum distance decoding. Define as := max{0,a;}. 

Theorem 1.1.2. Assume we use a code C for transmission over an operator channel. Let 
V E C he transmitted, and let 

u = nkiv)®E 

he received, where dim(£^) = t. Let p = (/(C) — k)+ denote the maximum number of erasures 
induced by the channel. If 

2{t + p)<d{C), 
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then a minimum distance decoder for C will produce the transmitted space V from the 
received space U. 



1.2 Rank-metric codes 

Let GF{q) be a base field and GF{q^) be an extension field of degree m. Consider an n- 
dimensional vector X = (xq, xi, G GF{q^)'^. The rank norm of X is defined by the 
maximum number of linear independent coordinates Xi over GF{q). The field GF{q'^) may be 
viewed as an m-dimensional vector space over GF{q). For any basis B of GF{q^) over GF{q), 
each coordinate of X can be expanded to an m-dimensional column vector over GF{q) with 
respect to B. Thus the rank norm of X is the rank of the m x n matrix over GF(q) obtained 
by expanding all the coordinates of X over GF{q). 

A rank distance between two vectors X, GF{q^)'^, denoted by dii{X,Y), is defined as 
the rank norm of the difference of X and Y , i.e. 

=rank(X-y). 

The minimum rank distance of a code over GF{q"^) is the minimum rank distance over 
all possible pairs of distinct codewords. The minimum rank distance of a code C of length n 
and cardinality M satisfies dji < du where dn is the minimum Hammimg distance of G. 
Therefore, by the Singleton bound we have: 

dR < dn < n — log^mM + 1. 

Codes that attain the equality are called maximum rank distance (MRD) codes. 
Now we consider the Gabidulin's construction of MRD codes [6]. 

Let hi € GF{q"^), i = 1, ...,n, n < m, be linearly independent over GF{q). Let d <n. Then 
the parity-check matrix 



H 



hi 

11 



ho 



h 



\hf-'^ h 



h 



[1] 

2 

,[2] 



[d~2] 

■2 



hn \ 



modr 



defines MRD code G with the rank distance d, where [i] = g' 
The generator matrix G of MRD code G with parity-check matrix H has the following form: 



G 



91 



92 

9? 



9n 
i/n 
i/n 



9n / 



where k = n — d + 1 and the elements gi,g2, ■■■,9n are linearly independent over GF{q). 
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Note that the size of such code is g'^C" d+i) ^ when m > n. Since rank(M) = rank(M-^) for 
each matrix M, we have that 



\C\ = q 



max{n,m} (min{n,m} — d+l) 



1.3 Organization of this work 

The rest of this work is organized as follows. In the next chapter we briefly overview the 
mathematical terminology which we use throughout this work. In Chapter 3 we consider the 
coding in projective space. First, we introduce the representation of subspaces in Vg{n) by 
constant-weight vectors and matrices in reduced echelon form. Then we provide an algorithms 
for computing the distance in Grassmannian and in Vq{n). Next we introduce a construction of 
constant-dimension codes based on optimal rank-metric codes and then, using this construction, 
we obtain a new codes in Vq{n). Finally using the connection between the Gaussian coefficient 
and the number of partitions we show the encoding methods for codes in Grassmannian. In 
Chapter 4 we discuss the future research goals. 
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Chapter 2 

Mathematical background 



In this chapter we introduce a few combinatorial notations and their properties which will be 
used throuhgout this work. 

Let GF{q) be a finite field of size q and V a vector space of dimension n over GF{q). 

Definition 2.1. Projective space (or projective geometry) of order n over GF{q), denoted 
by Vg{n), is a set of all subspaces of vector space V. 

Definition 2.2. Grassmannian, denoted by Gq{n, A;), is a set of all /c-dimensional subspaces 
of V, for any k < n. 

n 



Definition 2.3. Gaussian coefficient, denoted by 



k 



, is a q-analog for the binomial 



coefficient, also called a q-binomial coefficient or a Gaussian polynomial. A Gaussian coefficient 
is provided by 



n 
k 



k-l 

n 



1 



1 



J q i=0 

Theorem 2.4 [7]. The number of k- dimensional subspaces of an n dimensional vector space 
over GF{q) equals 



n 
k 



Lemma 2.5 |T]. The Gaussian coefficient 



n 
k 



satisfies 



-1 1 



n 
k 



< Aq 



k{n—k) 



Definition 2.6. A matrix is said to be in reduced echelon form if it has the following 
characteristics: 

1. All zero rows are at the bottom of the matrix 

2. The leading entry of each nonzero row after the first occurs to the right of the leading 
entry of the previous row. 

3. The leading entry in any nonzero row is 1. 

4. All entries in the column above and below the leading 1 are zero. 

Definition 2.7. A generator matrix for a fc- dimensional subspace U of GF{q)^ is a k x n 
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matrix over GF{q) of rank k whose rows span U . 

Note 2.8. Each fe-dimension subspace of GF{q)^ is uniquely defined by the generator matrix 
in reduced echelon form. 

Definition 2.9. A partition of a positive integer n is a representation of n as a sum of 
positive integers. Two sums which only differ in the order of their summands are considered to 
be the same partition. By convention, the summands of a partition are written from the largest 
to the smallest, for example, 10 = 3 + 2 + 2 + 2 + 1. The number of partitions of n is given by 
the partition function p{n). 

Example 2.10. p{5) = 7, the partitions of 5 are: 

1. 5 

2. 4 + 1 

3. 3 + 2 

4. 3 + 1 + 1 

5. 2 + 2 + 1 

6. 2 + 1 + 1 + 1 

7. 1 + 1 + 1 + 1 + 1 

Theorem 2.11. [S]. For all n > 0, the partition function p{n) is less then or equal to the 
(n + l)st Fibonacci number Fn+i- 

Definition 2.12. The representation of a partition by a diagram of dots called a Ferrers 
diagram. In the diagram, each term of partition represented by a row of dots, the terms are in 
descending order, with the largest at the top. 

Example 2.13. The partition (5,4,2,1) of 12 is represented by following Ferrers diagrams: 



Then the coefficient a; is the number of partitions of / whose Ferrers diagrams fit a box of 
size k hy n — k. 



Theorem 2.14 [7]. Let 
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Chapter 3 

Coding in projective space 



This chapter describes the main part of our ongoing research. We start by reviewing the back- 
ground and the related work. Next, we provide the novel representation of codewords in pro- 
jective space using binary vectors and matrices in reduced echelon form. We then use this 
representation to solve the following three problems of coding in projective space: comput- 
ing the distance between subspaces, constructing of constant-dimension codes, and encoding of 
constant-dimension codes. 



3.1 Background and related work 

Let V be an n-dimensional vector space over GF{q) and let Vq{n) be the projective space of 
order n over GF{q). It is known ^ that the dimension distance between any U,W G Vq{n) 
defined by 

d,{U, W) = dim(C/) + dim(VF) - 2dim([/ n W) (3.1) 

is a metric for the space Vq{n). A g-ary {n,M,d) code in projective space C is a subset of 
Vq{n) with size M and minimum distance d defined as 

d= min d^(X,Y). 

For a code C in Vq{n), the orthogonal complement code of C, C"*- is defined by C"*- = {X^ : X € C} . 
Etzion and Vardy in |9] proved that if C is an (n, M, d) code then C"*- is an (n, M, d) code. 

Let < n be any positive integer and let Gq{n, k) be the Grassmannian over GF{q). A g-ary 
(n, M, 25, k) constant- dimension code is a subset of ^g(n, k) with size M and minimum dimension 
distance 26. Note that it follows from (|3.1I) that the minimum distance of a constant-dimension 
code is always an even number and 1 < 5 <k. 

Let CC Qq{n,k) be an (n, M,26, k) constant-dimension code. Xia and Pu [10] proved that 
the orthogonal complement code of C, C"*- = {X-^ : X € C} is an {n, M,26,n — k) constant- 
dimension code. 
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3.1.1 Bounds on the size of the codes 



Denote by Aq{n,d) the maximum number of codewords in an {n,M,d) code in projective 
space, and denote Ag{n, 26, k) the maximum number of codewords in an (n, M, 26, k) constant- 
dimension code. An (n, Mi, d) code in projective space and an (n, M2, 26, k) constant-dimension 
code are said to be optimal, if Mi = Aq{n,d) and M2 = Aq{n,26,k), respectively. 

As was shown in pO], Aq{n,26,k) = Aq{n,26,n — k), therefore it is sufficient to determine 
Aq{n,26,k) for k < n/2. 

Let t = [^-^\ ■ Koetter and Kschischang [1] estabhshed the following sphere-packing and 
sphere-covering bounds: 



Aq{n,26,k) < 



\0q{n,k)\ 

\Sn,k{t)\ 



n 
k 





q 




' k ' 




n — k 














i=0 


i 


q 


i 


q 



(3.2) 



Aq{n,26,k) > 



\Gqin,k)\ 



n 
k 



\Sn,k{5-l)\ 5-1 


' k ' 




n — k 










i=0 


i 


g 


i 


g 



(3.3) 



where \Sn^k{t)\ is the volume of a sphere of radius t in Qq{n,k), which is well-defined since 
Qq{n,k) is distance regular. 

Koetter and Kschischang |1] also developed the Singleton type bound: 

Aq{n,26,k) < 



n-6 + l 
k-6 + 1 



(3.4) 



that is usually stronger than the sphere-packing bound (|3.2p . 
Xia in pT| showed a Graham-Sloane type lower bound: 



(9-1) 



Aq{n,26,k) > 



n 
k 



(g" - l)g"('^-2) 



but this bound is weaker than the bound 13.31 

In [T^ Wang, Xing and Safavi-Naini introduced the linear authentication codes. They showed 
that an [n, M, 26, k) constant-dimension code is exactly an [n, M, n — k,6] linear authentication 
code over GF{q). They also established an upper bound on linear authentication codes, which 
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is equivalent to the following bound on constant-dimension codes: 



Aq{n,26,k) < 



n 

k-5 + l 



k 

k-6 + 1 



(3.5) 



This bound was proved by using a different method by T.Etzion and A. Vardy in [9]. This 
method based on bounds on anticodes in Grassmanian. In [TO] was shown that the bound (|3.5p 
is always stronger then the Singleton bound (13.41) . Furthermore, it was proved that the codes 
known as Steiner structures attain the bound (|3.5|) . 

In [HI [To] the following Johnson- Type bounds were proved: 

1 



- 1 



Aq{n-l,2S,k - I), 
jAg{n-l,26,k). 

Using these bounds recursively, and combining with the observation that Aq{n,2d,k) 
for all k < 2(5, the following bound is obtained [9l [TO]: 



Ag{n,25,k) < 
Aq{n,25,k) < 



Aq{n,25,k) < 



g"-^ - 1 
r,k-l T 



q-n—k+S _ ^ 
- 1 



The upper and lower bounds on Aq{n, 25, k) when 5 = k were considered in [9]: 

- 1 



Aq{n,2k,k) < 



1 



1, if k\ n, 



Aq{n,2k,k) 



1 



if A; I n, 



qn _ gklqr _ 1) _ 1 

Aq{n, 2k, k) > , where n = r (mod k). (3.6) 

Etzion and Vardy [9j also obtained the following upper and lower bounds on Aq{n,d): 

Aq{n,d)> 



n d—1 j 

EEE 

A:=0i=0i=0 



n 


- k 




k 




n 


j 


— i 


1 


i 


q 


k 



Aq{n,2e + l)<r, 

where /* = max {Dq + Di + • • • subject to the following 2n -|- 2 linear constraints: 



^ c{i + j,i,e)Di+j < 



j=-e 



n 
i 
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and Di < Ag{n, 2e + 2, z), 

for all < i < n, where Di denote the number of codewords with dimension i and c{k, i, e) 
denote the size of the set {X : ds{X,Y) < e, dimX = i} for a /c-dimensional subspace Y. 



code in Gq{n, k) with M 



1 


' k ' 


t 


- 1 



and 6 = k — t + 1. Also they constructed a code 



3.1.2 Construction of codes 

Koetter and Kschischang in gave a construction and decoding procedures of Reed-Solomon like 
constant-dimension codes. They showed that these codes are nearly Singleton-bound-achieving. 

Etzion and Vardy [S] introduced a construction of codes in Gq{n, k) based on Steiner structure, 
that attain the bound (|3.5|) . They proved that any Steiner structure Sq{t, k, n) is an (n, M, 26, k) 

n 
t 

with 6 = k that attains the lower bound (13.6(1 . Etzion and Vardy [9] developed computational 
methods to search for the codes with a certain structure, such as cyclic codes, in Vg{n). 

Manganiello, Gorla and Rosenthal in [13] showed the construction of spread codes, i.e. codes 
that have the maximal possible distance in Grassmannian, and provided an efficient decoding 
algorithm for spread codes. 

Another constructions of codes in Grassmannian, given in [21 El 15) based on connection 
between rank-metric codes and constant-dimension codes. They will be discussed again in 
Section 3.4. 



3.2 Representation of codewords 



Let U G Vqiji) be any subspace of dimension k. Then there is a unique basis B = 
for [/, hi G GF{qy\ such that if we consider all the elements of i? as a rows of matrix U' in a 

/ h \ 

is a k X n matrix in reduced echelon form. Therefore, 



certain order, we get that U' 



U' is the generator matrix for U in reduced echelon form. In the sequel we will use a mixed 
language of subspaces and its generator matrices and denote by the same letter a subspace and 
its generator matrix in reduced echelon form. 

The following definition shows the connection between fc-dimensional subspaces in Vqin) and 
binary vectors of length n and weight k: 

Definition 3.2.1. Let U G 'Pq{n) be a subspace of dimension k. Consider U as a /c x n 
generator matrix for U in reduced echelon form. Define identifying vector Uid of length n and 
weight k such that if the leading entries of U are in columns ii,i2,---,ik G M, then for all 
1 < J < 

1, j = k,i2, ■■■,ik 
0, otherwise 



iU^d) 



This definition coincides with the definition of rank- and order-preserving map <p from Vq{n) 
onto the lattice of subsets of an n-set, given by D. Knuth in [H] and by S. Milne in [15], if we 
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consider an identifying vector Ui^ for some A;- dimensional subspace ?7 as a characteristic vector 
of a A;-subset. 

Note, that for each A;-dimensional subspace of GF(q)'^ exists only one identifying vector, but 
for each identifying vector there exist > 1 subspaces. 

Example 3.2.2. Let n = 4, A: = 2 and g = 2. Let U = span{(l, 0, 0, 0), (0, 1, 0, 0)} . Then, 
Uid = 1100. Now consider another /c-subspace = span {(1, 0, 0, 0), (0, 1, 0, 1)} . Obviously, 
U U^, but = Uid = 1100. Moreover, there are 2^ /c-subspaces that match this identifying 
vector. 

In the following we describe the construction of a family of subspaces corresponding to a 
binary vector. 

Definition 3.2.3. Let v G GF{2)'^ be a vector of weight k. Then echelon Ferrer form of v, 
denoted by EF{v), is the k xn matrix in echelon reduced form with leading entries in columns 
indexed by nonzero entries of v and " • " in the "arbitrary" entries. 

Example 3.2.4. Let v = 0110100. Then 

/oio.o..\ 
1 • • • . 

vooooi../ 

Lemma 3.2.5. Let v be a binary vector of length n and weight k with I's on places 
1 < ii, < n . Then, the number of /c-subspaces U over GF{q), such that Uid = 'u is equal 

to q~ '''V^^ +kn-{h+i2+---+ik) 

Proof. The number of such /c-spaces is equal to , where y is the number of dots in EF{y). 
The value of y is given by: 



EF{v) = 



{12 + 2(z3 -i2-l) + 3{u -i3-l) + ... + k{n - ik) 

-(1 + 2 + ... + fc - 1) - (ii + ... + ik) + kn 

k{k-l) , ... 
\- kn - [ii +i2 + ... + 



Corollary 3.2.6. Since the number of all A;-subspaces is equal to 



n 
k 



we have 



n 
k 



k(k-l) 



+kn 



l<ii<...<«fc<n 



(ii+.-.+ifc) 



3.3 Computation of the distance in Grassmannian 

In this section we present the algorithm for computing the distance between subspaces in Grass- 
mannian. We first present a few lemmas that will be used in this algorithm. We consider the 
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case q = 2, but as we will note at the end of the section, all the results can be easily generalized 
to each q. 

Let U and W two subspaces in Q2{'n, k). Recall that the distance between U and W is defined 

by 

d,{U, W) = 2k- 2dim(i7 n W). 



Now denote hy U *W the concatenation of U and W, i.e.: 

U*W - 



U 
W 



Lemma 3.3.1. 

dim(C/ nW) = 2k- rank{U * W). 

Proof. rank(;7 *W) = rank(C/) + rank(W^) - dim(i7 nW) = 2k- dim(;7 n W). 
U. 

Lemma 3.3.2. 

d,{U, W) = 2rank(?7 *W)- 2k, 

where W and U are also considered as the generator matrices in reduced echelon form for 
W and U, respectively. 
Proof. By Lemma 3.3.1, 

d,{U, W) = 2k- 2dim(f/ nW) = 2k-Ak + 2rank(C/ *W) = 2rank(C/ *W)- 2k. 

U. 

Corollary 3.3.3. 

d,{U, W) = 2d^ rank{U *W) = k + d. 

U. 

Let vi,V2 be two binary vectors of the same length. Denote by dH{vi,V2) the Hamming 
distance between vi and V2, i.e. the number of different entries of these two vectors. 

Let U and W be two /c-subspaces with dHiUid,Wid) = 2s. Define three sets /, J, L in the 
following way: 

/ = (1 < < n| [Uid]i^ = 1 and [Wial^^ = 0} , 

J={l<jt<n\ [W^d]j, = 1 and [Uid]j^ = o} , 

L={l<lt<n\ [Uidli^ = 1 and [Wid]i^ = l} . 

Suppose w.l.o.g. that ii < ji. We say that a row v in A E {I,J,L} if its leading entry in a 
column indexed by A. Note that |/| = | J| = s and \L\ = k — s. 



15 



Let Ui be the s x n sub-matrix of U that consists of the rows in / and let Wj be the s x n 
sub-matrix of W that consists of the rows in J. Denote by [U * Wj] and [W * Uj] the reduced 
echelon form of U * Wj and W *Ui, respectively. Let U and W he {k — s) x n submatrices of 
[U * Wj] and [W * Uj], respectively, with rows in L. 

Lemma 3.3.4. 

vank{U *W) = k + ^dHiUid, Wid) + rank(t/ + W), 

where U and M are defined above. 
Proof. 

Consider the process of bringing to reduced echelon form oi U *W . We have the following 
stages: 

U 



U 

w 



Wj 
W\Wj 



' [U * Wj] ' 




' [U * Wj] ' 


w 




w + u 



Therefore, rank{U*W) < raxik{[U *Wj]) + rank{W + U). Note that all the rows of [U*Wj] 
are in / U J U L while by definition of W and U, all the rows of W + U are in [n] \ I U J U L. 
Thus, W + U ^ span{[i7 * Wj]} , and we have the following equivalence: rank(C/ * W) = 
rank([?7 * Wj]) + rank(l^^ + U). Since rank([?7 * Wj]) = \I\ + \L\ + \J\ = k + ^niUid, Wid), we 
have that rank(C/ *W) = k + ^dniUid, Wid) + rank{W + U). 

U. 

Corollary 3.3.5. If d,{U, W) = 2d and dniUid, Wid) = 2s, then d>s. 
Proof. By Lemma 3.3.2 and Lemma 3.3.4: 

2d = ds{U, W) = 2 rank(?7 * 1^) - 2A; > 2A; + 2s - 2A; = 2s. 



Lemma 3.3.6. Let U and W be two A;-subspaces with dH{Uid,Wid) 
Then 

ds{U, W) = 2d^ rank(C/ + W) = d-s, 

where s < d. 

Proof. By Corollary 3.3.3 and Lemma 3.3.4: 



2s. 



ds{U,W) = 2d 



k + d = rank(;7 *W) = k + s + rank(?7 + W) 
rank{U + W)=d-s. 



Corollary 3.3.7. Let U and W be two /c-subspaces with dniUid, W^d) = 0. 
Then 

ds{U, W)=2d^ rank(?7 + W) = d. 
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Corollary 3.3.8. 

ds{U, W) = dniUu, Wid) + 2rank(t/ + W). 

U. 

Note that it is easy to generalize all the above for all q. The only difference is that instead 
of addition we have to use subtraction, or in other words, instead of a + 6 we need write 
a + {q — 1)6, where a,b G GF{q). That is, in each place where appears rank(C/ + W) we put 
now rank(?7 + (g — 1)1^) • 

3.3.1 Algorithm for computing the distance between two k - subspaces. 

Using the Corollary 3.3.8. we can easily determine the distance between any two subspaces 
having the same dimension. 

Let U and W be two A;-subspaces over GF{2) with the generator matrices in reduced echelon 
form U and W , and identifying vectors Uid and Wid, respectively. 

1. Calculate the Hamming distance between Um and Wi^, and denote 

dH{Uid,Wu) = 2s, 

since the distance is even. 

2. Define the following two sets: 

I={it\[Uid]i, = l^^d [Wid]i,=Q], 

J = {jt\ [Wid]j^ = 1 and [Uid]j, = O} . 

Suppose w.l.o.g. that ii < ji (otherwise we swap the roles of U and W). Note that 
|/| = \ J\ = s. Add each row in J with the leading entry jt to the rows of U with 1 in the 
jt coordinate, and add each row in I with the leading entry it to the rows of W with 1 in 
the it coordinate, in decreasing order of the elements of / U J. Denote by U and W the 
resulting k — s x n matrices without rows in / U J. 

3. Calculate the rank olU + W. 

4. Output dsiU, W) = 2s + 2rank(L/ + W). 

Now consider complexity of the algorithm above: 

• comparing two identifying vectors : n 

• obtaining the reduced echelon form < sk{n — k) + s{k — s){n — k) = {2ks — s^){n — k) 

• computing the sum {U + W): note that we only need to consider submatrices of size 
{k — s) X [n — k — s) that correspond to common zero entries in identifying vectors: 
{k — s){n — k — s) 
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• computation of rank oi {k — s) X (n — k — s) matrix : 

min {(/c — s)^(n — k — s), {k — s){n — k — s)^} 

Thus, the complexity of above algorithm (assuming that k < n — k): 

(n + {2ks - s^)(n - k) + {k - s){n - k - s) + {k - sf{n -k- s)). 

Note that this algorithm may be useful for decoding in Grassmannian. That is, if given any 
/c- dimensional subspace W we want to find a codeword C with ds{C, W) < d, we have to consider 
only the codewords which satisfy the condition dH{W,C) < d, by Lemma 3.3.5, and then given 
a list of such codewords we need to apply the algorithm above and check if ds{W,C) < d for 
each codeword in the list. 

3.4 Computation of the distance in V2{ti). 

In this section we present the algorithm for computing the distance between subspaces in 7^2 (?^), 
which is the generalization of the results obtained in the previous section in G2in, k). As before, 
we start with a few lemmas. 

Let U and W be two k\- and A;2— subspaces, respectively, in V2{n) 

Recall that the distance between U and W is defined by 

d,{U, W) = ki + k2- 2dim(?7 n W). 

Lemma 3.4.1. 

dim(?7 ^W) = ki + k2- rank(?7 * W). 
Proof. rank(C/ *W) = rank(C/) + rank(VF) - dim(C/ r\W) = ki + k2 - dim(C/ n W). 

m. 

Lemma 3.4.2. 

d,(U, W) = 2rank(C/ *W)-ki-k2, 

where W and U are also considered as the generator matrices in reduced echelon form for 
W and U, respectively. 
Proof. By Lemma 3.4.1, 



ds{U,W) = ki + k2-2dim{UnW) 

= ki + k2- 2{ki + k2) + 2rank(?7 *W) = 2rank(?7 *W)- {ki + k2). 

U. 

Corollary 3.4.3. 
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d,{U, W)=d^ rank(i7 * W) 



ki + k2 + d 



Lemma 3.4.4. Let Uid and Wid be the identifying vectors for U and W, respectively. Then 

dniUid, Wid) = ki + k2- 2Uid ■ Wid, 
where "•" is a dot product. 

Proof. Since weight (?7i(i) = ki, weight (W^^d) = k2, and Uid • Wid is equal to the number of 
common nonzero entries of Uid and W^d, we have that dniUid, Wid) = ki + k2 — 2Uid ■ W^d- 
U. 

Lemma 3.4.5. ds{U,W) = dH {Uid, Wid) (mod 2). 

Proof. From the definition of the distance dg, it follows that ds{U,W) = {ki + A;2)(mod2). 
From the Lemma 3.4.5. dH{Uid,Wid){ki+k2) = (A;i+A;2)(mod2), hence ds{U,W) = dH{Uid,Wid){mi 



Let U and W he ki- and A;2— subspaces, with dH{Uid,Wid) = s. Let define three sets /, J 
and L in the following way: 



and 



/ = {1 < it < n| [Uid]i^ = 1 and [Wid]i^ = O} , 
J={l<jt<n\ [W,d],, = 1 and [Uid]j^ = o] , 

L={l<lt<n\ [Uid]i^ = 1 and [Wid]i^ = l} • 



Suppose w.l.o.g. that ii < ji. We say that row v in A E {I,J,L} if its leading entry is in a 
column indexed by A. Note that |/| + | J| = s and \L\ = Uid " Wid = (^i±|2l:i£. Denote by si the 
size of /, i.e. |/| = si, thus \L\ = ki — si. 

Let Ui be the si x n sub-matrix of U that consists of the rows in / and let Wj be the 
(s — si) X n sub-matrix of W that consists of the rows in J. Denote by \U * Wj] and \W * Ui] 
the reduced echelon form oiU * Wj and W * f//, respectively. Let U and W he {ki — si) x n 
submatrices of [U * Wj] and [W * Uj], respectively, with rows in L. 

Lemma 3.4.6. Let U and W he ki- and A;2— subspaces. Then 

rank(C/ * W) = h±h. + ^dniUid, Wid) + rank(t/ + W), 

where U and M are defined above. 
Proof. 

Consider the process of bringing to reduced echelon form of U *W . We have the following 
stages: 

U 



U 
W 



Wj 

W\Wj 



' [U * Wj] ' 




' [U * Wj] ' 


w 




w + u 
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Therefore, rank(?7*T^) < rank([?7* PFj]) + rank(PF + i7). Note that all the rows of [?7*PFj] 
are in /U JUL while by definition oiW and U, all the rows oiW + U are in [n]\/U JUL. Thus, 
W + U ^ span{[?7 * VFj]}, and we have the following equivalence: rank(C/ * W) = rank([C/ * 
Wj]) + iank{W + U). Sincerank([i7*W^j]) = |/| + |L| + |J| = dH{Uid,Wid) + ^^^^^^^^^f^^^^ = 
dH(U^d,w^^d)+kl+k2 ^ follows that rank(?7 *W) = + \dH{Uid, Wio) + rank(PF + U). 

m. 

Corollary 3.4.7. If ds{U,W) = d and dH{Uid,Wid) = s, then d>s. 
Proof. By Lemma 3.4.2 and Lemma 3.4.6: 

d = ds{U, W) = 2 rank(i7 *W) - ki - k2 > h + k2 + s - h - k2 = s. 



Lemma 3.4.8. Let U and W be two /c-subspaces with dH{Cid,Wi(i) = s. 
Then 

d,{U, W) = d^ rank(f7 + W) = 

where s < d. 

Proof. By Corollary 3.4.3 and Lemma 3.4.6: 



ds{U,W) = d 



{ki +k2) + d {ki + k2) + s ~ ~ 

^ = rank(i7 * W) = ^ h rank(?7 + W) 

rank{U + W) = 



Corollary 3.4.9. 

ds{U, W) = dniUu, Wid) + 2rank(t/ + W). 

U. 

Note that it is easy to generalize all the above for all q. The only difference is that instead 
of addition we have to use subtraction, or in other words, instead of a + 6 we need write 
a + {q — 1)6, for a, 6 G GF{q). That is, in each place where appears rank(J7 + W) we put now 
rank{U +{q-l)W). 

3.4.1 Algorithm for computing the distance between two subspaces in V2{n). 

Using the Corollary 3.4.9. we can calculate the distance between any two subspaces in V2(n). 

Let U and W be two ki- and A;2-subspaces with the generator matrices in reduced echelon 
form U and W , and identifying vectors Uid and Wid, respectively. 

1. Calculate the Hamming distance between Uid and Wid, and denote 

dH{Uid,Wid) = s. 
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2. 



Define the following two sets: 



J = 



[jA [Wid]j, = 1 and [Uid]j, = 



} 



Suppose w.l.o.g. that ii < ji (otherwise we swap the roles of U and W) and denote by si 
and S2 the size of / and J, respectively (note that si + S2 = s). Add each row in J with 
the leading entry jt to the rows of U with 1 in the jt coordinate, and add each row in / 
with the leading entry it to the rows of W with 1 in the it coordinate, in decreasing order 
of the elements of /U J. Denote by U and W the resulting {ki — si) x n matrices without 
rows in / U J. 

3. Calculate the rank oiU + W. 

4. Output ds{U, W) = s + 2rank(t/ + W). 

Now consider complexity of the algorithm above: 

• comparing two identifying vectors : n 

• obtaining the reduced echelon form < (siA;2 + ^2(^1 — si)){n — ki — S2) ■ 

• computing the sum {U + W): note that we only need to consider submatrices of size 
{ki — si) X [n — ki — S2) that correspond to common zero entries in identifying vectors: 
{ki - si){n -ki- S2) 

• computation of rank of (A:i — si)(n — fci — S2) matrix : 



Thus, the complexity of above algorithm (assume that the number of common ones is less then 
the number of common zeros, i.e. {ki — si) < {n — ki — S2)): 



Note that this algorithm may be useful for decoding in V2{n). That is, if given any subspace 
W we want to find a codeword C with ds{C,W) < d, we have to consider only the codewords 
which satisfy the condition dH{W,C) < d, by Lemma 3.4.7, and then given a list of such 
codewords we need to apply the algorithm above and check if ds{W,C) < d for each codeword 
in the list. 



min {(A;i - sif{n - ki - S2), (fci — Si)(n — ki — 82)^} 



{n + {n-ki- S2) [{sik2 + S2{ki - si)) + {ki - si) + {ki - si)^] . 
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3.5 Construction of codes in Pg(n) from rank-metric codes 



3.5.1 Codes for the random network coding channel based on rank- metric 
codes 

Codes in the rank-metric were studied by a number of authors [6l dU [171 Ell US] and have been 
recently proposed for error control in random network coding [ll[3l[2]. 

In this section we describe the construction of codes for the random network coding based 
on rank-metric codes, called the lifting construction, given by D. Silva, F. Kschischang and R. 
Koetter [2]. 

Definition 3.5.1.1. The lifting of a matrix x £ GF{q)^^™' is defined to be matrix X = [I x] 
obtained by prepending to x an identity matrix. Similarly, the lifting of rank-metric code 
C C GFiq^"^ is a code 

l[C] = {[Ix], xeC} 
obtained by lifting every codeword x £ C. 

Given any rank-metric code C C GF{q)^^"^ with minimum distance da it can be constructed 
a subspace code X = I[C] that is a constant-dimension code with minimum distance 2dR. 

The next proposition shows that this construction is "good" compared to the whole class of 
constant-dimension codes: 

Proposition 3.5.1.2. Let CC GF{q)"-^"^ be an MRD code with minimum distance dR. 
For any subspace code Q with minimum distance d{Q) > 2dR, we have < 4|C|. 



3.5.2 Definitions 

Our construction is a generalization of the lifting construction, based on presentation of constant- 
dimension spaces by constant-weight binary codes. 

Recall that the echelon Ferrer form oiv £ GF(2)", denoted by EF{v), is the kxn matrix in 
reduced echelon form with leading entries in columns indexed by nonzero entries of v and " • " 
in the "arbitrary" entries. 

Example 3.5.2.1. Let v = 0110100. Then 



EF{v) 



/ 1 • 
1.0 

V 1 



Definition 3.5.2.2. Let v £ GF(2)" and M G ^^^^^ 1 < < n - k is 

the index of the first one of the v. We say that M € EF{v) if the nonzero entries of M can be 
fit the point part of EF{v). 



Then 



Ml G EF{v) and Ma ^ EF{v). 

Definition 3.5.2.4. Let v £ GF{2Y and M £ EF{v). Define the matrix EF{v[M]) be 
the matrix in reduced echelon form that resulted by putting the matrix M in the point part of 





( 


1 


1 






( 


1 


1 




Example 3.5.2.3. Let v = 0110100, Mi = 




1 


1 









1 


1 






\ 





1 






K 


1 


1 
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EF{v). 

Example 3.5.2.5. Let v = 0110100, Mi 



/I 1 1 

1 1 I from the previous example. Then 
V 1 1 



EF{v[Mi]) 



/ 1 1 1 1 
110 10 

V 1 1 1 



Definition 3.5.2.6. Let v G GF{2y and S = {Mi, Ma, ...,MJ, such that Mi G EF{v) for 
all 1 < i < t. Define the set of matrices EF{v[S]) in the following way: 



EF{v[S]) = {EF{v[Aml,. 



3.5.3 Construction in Grassmannian 

Let be a constant-weight code of length n, constant weight k and minimum Hamming distance 
dn = 25, that includes the vector 111. .. 1000. ..0. 

For each w G W we construct the echelon Ferrer form, EF{w). Let be an optimal rank- 
metric code with the minimum distance dR = 5, such that all its codewords in EF{w). For each 
w we construct a subset C{w) of code C as follows: 

C{w) = EF{w[C^]). 

Now define code C by 

c= u CH. 

Lemma 3.5.3.1. The constructed code C is constant-dimension code. 

Proof. The generator matrix for each codeword in C is a matrix in reduced echelon form 
with k nonzero rows, therefore all the codewords have the same dimension k. 

Lemma 3.5.3.2. The minimum distance of the constructed code C is ds = 25. 

Proof. By Corollary 3.3.5, for each V,U eC, ds{V, U) > duiVid, Uid)- Therefore, ds > dn = 
25. Now let V,U e C{w), for some w eW, such that dR{V,U) = rank(F -U) = 5. By the 
Corollary 3.3.7, since duiVid, Uid) = 0) it follows that ds{V, U) = 25. Therefore, ds = 25. 

By using the construction, we provide an alternative proof of the following theorem which 
first appeared in [H]. 

Theorem 3.5.3.3. Let n = r( mod k). Then, for all q, we have 

A,{n,2k,k) > ^ '^[^ ' . 

g*= — 1 

Proof. Let n = sk + r, where < r < fc. Let C = {ci,C2, ...,Cs} be a constant-weight code 
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with the minimum distance 2k with the codewords of the following form: 

ci = 11.. .10000000000000000000.. .00 

k n—k 

C2 = 00. ■ m 1 ■ ■ . 1 00000000000000. ■ .0 o 

k k n—2k 




{s-2)k k n-{s-l)k 




(s-l)fe k r 



For each 1 < i < s - 1 let Q G GF{qf''-'^ be the Gabiduhn MRD code with minimum 
distance k. Note that each codeword of Q is in EF{ci) and \Ci\ = qi'^-i-k){k~k+i) _ ^n-tk p^j. 

s 

Cs we take Cg that include only k X r zero matrix, thus \Cs\ = 1- Let C = U £^F(cj[Cj]). By 

i=l 

construction, this is a constant-dimension code with minimum distance 2k. Now we calculate 
its size: 



^|Q|=^g"-^^ + l 



=1 i=l 
s-1 



i=l 



hs-\-h _|_ qh q: 



.k _ I 



n r/k+'f' 



,k-k{s-l) _|_ 



q^ - 1 



q^ - 1 



Theorem 3.5.3.4. Let C C GF[q)™'^^^ be rank-metric code with minimum distance (Ir 
such that all its codewords of the following form: 



XXX 



\ 



X 





X X 

X 



X X 
X X 



... X ... X 
VO... X ... X J 



where x is any element of GF{q), and denote by Ij the number of zeros at the beginning of 
the (m — j)ih. row, < j < m — 1, and denote by ji the number of zeros at the end of the zth 
column, 1 <i <n. 

Then the size of the code is 



\C\<q 



-d+l n-d+1 
3 = 1 i=l 
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Note that the exponent of the above expression is the minimum between the number of 
nonzero entries without nonzero entries in the first d — 1 rows and number of nonzero entries 
without nonzero entries in the last d — 1 columns. 

Proof. Suppose that there is a code C over GF{q) of the above form with dim(C) = t + 1, 

{m-d+l n-d+l \ 

^ (n — lj), ^ {m — ji) >. Let B = {bi,b2, ■■■,bt+i} be a basis for C. 
i=i i=l J 

Then there is a nontrivial linear combination Yli=i ^nh = M of the following form: 

fx ... X X X ... X \ 



X . 


. X 




X 


X 


. 


. 


.. 


. 






... ... 
VO... ... 0/ 



m—d+l 

a t = {n — Ij), or 



/O ... X ... x\ 



. 


. 


X . 


. X 


. 


. 


X . 


. X 



... X ... X 
\0 ... X ... X j 



n-d+l 

if i = (i^ ~ ^-nd in any case this is a matrix of rank d — 1, that is contradiction, 

i=l 

since minimum rank-distance is d. 
U. 

Conjecture 3.5.3.5. There is a rank-metric code C G GF{q)'^^^ with minimum distance 
dR of the above form, that achieve the bound of Theorem 3.5.3.4. 

Note that the largest part of the resulting code is obtained from the constant- weight codeword 
of the form 111. ..1000. ..00, and its size is ^ that is equal to the size of the MRD 

rank-metric code in GF{q)^^^^~^\ The sizes of other parts of the resulting code are bounded 
by the g("~'^)('^~'^+^)~''^each. Therefore the contributions of these other parts becomes more 
significant for the lower values of 5. Thus in the following we will focus on the case 5 = 2. 

For 5 = dn = 2, the following theorem shows that there is a code that achieve the bound 
of the Theorem 3.5.3.4, therefore for this minimum rank distance we have that the Conjecture 
3.5.3.5 is true, but before this we rewrite the bound of Theorem 3.5.3.4 in another way: 
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mn— ^ ij— max{m,n} 

\C\ < q 

where the number of zeroes at the beginning of the last row is ii, in the one before the last 
is i2, and in the first row with zeroes is ig, and s <m and ij < n, for all 1 < j < s. 

Note that for the case s = m and ij < n we can consider the code C C GF{q)"^^'^ , where 



n = n 



and 



m 



= m-\{ij -.ij =n}\. 



Theorem 3.5.3.6. Let C G GF{q)'^^'^, m < n, he the Gabidulin rank-metric code with 
minimum distance (Ir = 2. Then there is the subcode C of C such that each codeword in C' of 
the form 

/ X ... X X X ... X \ 



tXJ ... tly it/ ... it. it. 

... X ... X X 



... X ... X 
V ... X ... X J 

where x is any element of GF{q), and the number of zeros at the beginning of the last row 
is ii, in the one before the last is 12, and in the first row with zeros is is, s < m. Then the size 
of the code is 



\C'\ = q . 

Proof. Note, that since m < n, we swap the roles of n and m in the Gabidulin's construction, 
and transpose the codewords of Gabidulin code. The parity-check matrix H of the code C is of 
the following form: 

H = hi h2 ... hjn^ , 

where hi G GF{q^), i = 1, ...,m, m < n, are linearly independent over GF{q). Let a be a 
primitive element of GF{q'^). The field GF{q^) may be viewed as an n-dimensional vector space 
over GF{q) with basis B = |l,a,a^, ...,q;"~^} . We choose hi = a*"^, 1 < i < m. Then 



H=il 



a a 



a 



TO — 1 



Now we can construct the (m — 1) x m generator matrix of the following form: 
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/ a 1 ... \ 
q; 1 ... 



G = 



... a 10 
y ... a I J 

Note that each codeword can be represented as a hnear combination of the rows of G over 
GF{q'^). Therefore, the general element c G C is of the form: 

/ aia \ 

ai + a2cx 
02 + asa 

a-m-2 + am-lOi 

where a-i € GF{q^^), 1 < i < m— 1. Now we expand each coordinate of c to an n-dimensional 
row vector witli respect to the basis B and obtain the m, x n matrix. If c G C", then we have 
the following system of equations, where the variables are coefficients Oj, 1 < i < m — 1: 



Qm-i = OO...Oxx...x 

ii 

am-2 + CLm-iOi = OO...Oa;x...x 

12 

< am-s + am-s+ia = OO^a;a;...x 

is 



aia = xx^^ 

n 

where x is any element of GF{q). The size of the code C is equal to number of solution of 
the system above. 

There are g""'! possibilities for a^_i, g"'~*2 possibilities for am-2, ■■■ ,q"'~^" possibilities for 
am-s and possibilities for all Oj, l<i<m — s — 1. Therefore, 



\C'\ = . . . qn-is^^qn^^im-s-l) 



ns—'^ij+nm—ns—n nm—'^ij—n 
q ■'■=1 = q 



Prom the Theorem 3.5.3.6 and Lemma 3.2.5. we obtain the following corollary: 
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Corollary 3.5.3.7. Let f be a binary vector of length n and weight k with I's with the 
indices 1 < ii,...,ik < n. Then, the size of optimal rank-metric code with the minimum 
distance (Ir = 2 such that all its codewords in EF{v) is equal to g^"-'' \ -(«i+J2+-+*fe)-6^ 
where h = max {n — /c — ii + 1, fc — (# of ones before the last zero)} . 

■ . 

Example 3.5.3.8. Let g = 2, n = 5, = 2 and 4 = 4. Let C = {ci = 11000, C2 = 00110} 
be the constant-weight code. Then by Corollary 3.5.3.7, rank-metric code Ci for ci is of size 
8, and the rank-metric code C2 for C2 is of size 1. Thus we got the constant-dimension code of 
length 5, dimension 2, minimum distance 4 over GF{2) of size 9, (5, 9, 4, 2) code, that is by 
[9j is an optimal code. 

Example 3.5.3.9. Let g = 2, n = 6, fc = 3 and ds = 4. We summarize the results for these 
parameters in the following table: 



i 


Ci- codeword of constant-weight code 


size of corresponding rank-metric code Cj 


1 


111000 


64 


2 


100110 


4 


3 


010101 


2 


4 


001011 


1 



Therefore we constructed the constant-dimension code of length 6, dimension 3, minimum 
distance 4 over GF{2) of size 71, thus ^2(6,4,3) > 71, which improves the bound in |9]. 

Example 3.5.3.10. Let g = 2, n = 7, A;i = 4, ^2 = 3 and ds = 4. We summarize the results 
for two constructions, ki- and /c2 -dimensional codes, respectively, in the following table: 



i 


Q- codeword of 


Cj- codeword of 


size of corresponding 




constant-weight code 


constant-weight code 


rank-metric code Ci 




with weight 3 


with weight 4 




1 


1110000 


1111000 


256 


2 


1001100 


1100110 


16 


3 


0101010 


1010101 


8 


4 


0010110 


1001011 


2 


5 


0100101 


0101101 


2 


6 


0011001 


0110011 


4 


7 


1000011 


0011110 


1 



Therefore we have that ^2(7,4,3) = ^2(7,4,4) > 289. 



Example 3.5.3.11. Let g = 2, n = 8, fc = 4 and ds = 4. We summarize the results for 
these parameters in the following table: 
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i 


Cj-codeword of 


size of corresponding 




constant-weight code 


rank-metric code Q 


1 


11110000 


4096 


2 


11001100 


256 


3 


10101010 


64 


4 


10010110 


16 


5 


01011010 


16 


6 


01100110 


16 


7 


00111100 


16 


8 


11000011 


16 


9 


10100101 


16 


10 


10011001 


16 


11 


01010101 


8 


12 


01101001 


32 


13 


00110011 


4 


14 


00001111 


1 



Therefore we have that ^2(8,4,4) > 4573. 



Open problem 3.5.3.12. What type of a constant-weight code results in optimal constant- 
dimension code constructed as above? 

Note, that the maximal code with n = 8, A; = 4 and = 4 of the last example obtained 
from all codewords of the weight 4 of extended Hamming code with the following parity-check 
matrix: 

/ 1 1 1 1 1 1 1 1 \ 
1 1 1 1 , , 
1 1 1 1 

^ 1 1 1 1 y 

Conjecture 3.5.3.13. An optimal constant-dimension code with n = 2'' , A; = 4 and c?s = 4 
obtained from all codewords of weight 4 of extended Hamming code with the parity-check matrix 
of the similar to (*) form. 
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3.5.3.1 Construction of rank-metric codes of a given form with minimum rank 
distance 2 

Let X he n X m matrix over GF{q) of the following form: 

fx ... X X X ... X \ 



X 



X 





X X 

X 



X 
X 



X 
X 




V 








X 
X 



X 

X J 



where x is any element of GF{q) and the number of zeros at the beginning of each row does 
not decrease from the top to bottom. We want to construct directly a rank-metric code with all 
codewords in the same form. So far we solved this problem only for a relatively small subset of 
parameters. We summarize the results in the following examples: 

Example 3.5.3.1.1. Let n = m = 3. We have that each codeword in the code C = 
span {ci, C2, C3, C4, C5}, with minimum distance 2, where 




C2 



/ 1 
10 
V 



C3 






( ° 





1 \ 




( ° 





1 \ 


C4 = 


1 





















I 





0) 




I 


1 


0/ 



of the form 

^ X X X 

x X x \ (**). 
\ X X 

Moreover, each form that have more zeros than (**) can be obtained by choosing suitable 
elements of basis above. 

Example 3.5.3.1.2. Let m = 3, n = 4. Than such basis with given parameters is: 



ci 



/ 1 1 \ 
1 

V y 



,C2 



/ 1 \ 
10 

V / 



,C3 



/ 1 \ 


V 1 / 



,C4 



/ 1 \ 
10 

V / 



10 
C5 = I 
10 



,C6 





30 



and from these matrices can be constructed a matrix of the following form: 

^ X X X X ^ 

tAj tXj J 

\ Q X X X J 

and therefore, each matrix with more zeros. 

Example 3.5.3.1.3. For m = 3 and general n, the following basis can be used to constructed 
all matrices with zero entry in the left bottom: 



/o 





. 


. 





1 







. 


. 


1 





\ 







. 


. 


1 





o\ 








. 


. 








;)■(; 





. 


. 





1 










. 


. 











Vo 





. 


. 













. 


. 








J 







. 


. 








1 / 









1 \ 
10 
/ 






. 


. 1 














. 


. 














. 


. 





1 






10 
10 










1 





. 


. 














. 


. 














1 . 


. 












. 


. 










^0 





. 


. 








1 \ 





. 


. 
















1 . 


. 














. 


. 










^0 


1 


. 


. 











Example 3.5.3.1.4. For m = n = 4, the basis is: 



ci 



/ 1 1 \ 

1 



y y 



, C2 








1 








\ 




/ 





1 








\ 




/ 


1 





















1 
































1 













































































1 






















v 














/ 




\ 














/ 




I 














) 



C5 = 



/ 


1 















/ 


1 















( 











1 


\ 




/ 











1 


\ 








































1 





























































































1 























































V 














) 




I 











1 


) 




\ 














/ 




I 








1 





/ 



eg = 



/ 











1 






/ 








1 









/ 














\ 










1 














1 























1 
















, ClO = 














, Cll = 



















1 












1 





















1 







V 














) 




I 














) 




I 





1 








/ 
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3.5.4 Construction of codes in the projective space from the constant-dimension 
codes based on rank-metric codes 

The general construction is as follows. Let C be a constant-dimension code of length n, dimension 
k and minimum distance 26. Let Q be some (n — 1)- dimensional subspace. Define two subcodes 
Ci, C2 C C by Ci = {c G C| c C Q} and C2 = {c G C : v e c} D Q, for any vector v ^ Q of 
length n. Then Ci and C2 are constant-dimension codes of length n — 1, minimum distance 25, 
dim(Ci) = k and dim(C2) = k — 1. Note that the distance between ci G Ci and C2 G C2 is 
ds(ci, C2) > A; + — 1 — 2{k — 5) = 25 — 1. Now we define new code C = Ci U C2. C is of length 
n — 1, minimum distance 26 — 1 and size |C[ = \Ci \ + IC2I 

Let C be a constant-dimension code constructed from constant-weight code and rank-metric 
codes. Let Q be an n — 1-dimensional subspace with the following generator matrix: 

/ 1 ... \ 
1 ... 



... 1 

... 1 y 

n 

and let v = 100.. .001 or 0^100.. .01. Then C = {c G C : c C Q} u ({c G C : v G c} n Q) is 

k-l 

an (n — 1, |C|,2(5 — 1) code in projective space. 

Example 3.5.4.1. Let C be the (6, 71, 4, 3) constant-dimension code from Example 3.5.3.9. 
Let consider Ci. Note that c G Ci iff c G C and the last column of generator matrix of c is zero. 
Therefore Ci is obtained from two constant-weight codewords: ci = 111000 and C2 = 100110. 
The MRD code associated with ci is set of 3 x 3 matrices with rank distance 2, thus its subcode 
with zero last column is subspace subcode, hence by [16] this is a MRD code of size 8. The 
second codeword, C2 gives only one codeword in Ci. Thus we get that |Ci[ = 9, that is an 
optimal (5, 9, 4, 3) code. Now consider C2 = {c G C : 001001 G c} PI Q. The code C2 is obtained 
from two constant-weight codewords: ci = 111000 and C4 = 001011. Let be rank-metric 
code associated with ci, consider its subspace subcode with zero last row in its generator matrix 
in reduced echelon form. By [16], this is a MRD code of size 8. Now we sum all codewords of 
this subspace subcode with any codeword of Cci with the last row of the form 001. Thus we 
get 8 codewords of C2. The second codeword C4 gives one codeword in C2, thus IC2I =9 , and 
it is an optimal (5, 9, 4, 2) code. 

As a result we get code Ci U C2 of minimum distance 3 of size 18, thus ^2(5, 3) > 18. 

Example 3.5.4.2. Let C be the (8, 4573, 4, 4) constant-dimension code from Example 
3.5.3.11. Note that Ci is obtained from ci,...,C7 and equal to the constant-dimension code 
(7, 289, 4, 4), since the MRD rank-metric code associated with 11110000 is the set of square 4x4 
matrices and its subspace subcode is also MRD code [16]. Thus, |Ci| = 289. Now consider C2 = 
{c G C : 10000001 G c} n Q. The code C2 is obtained from codewords ci, C2, C3, C4, cs, cg, cio- 
Let Cci be rank-metric code associated with ci, consider its subspace subcode with zero first 
row in its generator matrix in reduced echelon form. By [16], this is a MRD code of size 256. 
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Now we sum all codewords of this subspace subcode with any codeword of Cci with the first 
row of the form 0001. Thus we get 256 codewords of C2. We use the construction of basis of a 
given form from example 3.5.3.1.4 in order to calculate the size of rank-metric code for each Cj, 
i G {2,3,4,8,9, 10}. The results are summarized in the following table: 



i 


Cj-codeword 
of 

constant- weight 
code 


generator matrices of corresponding 
rank-metric code (the coefficients are 
in {0,1}) 


size of 
corresponding 
rank-metric code 


2 


11001100 


Cy -|- a(co + C'^) -|- b(cA + Ct;) -|- dfcc; -|- Cp.) 

or 

C8 + a(c2 + C3) -1- 6(c4 -1- C5) -1- (i(c5 -1- ce) 


16 








3 


10101010 


(c5 + ce) -f- a(c2 -1- C3) 5c9 


4 


4 


10010110 


(c5 + ce) 


1 


8 


11000011 





1 


9 


10100101 


a(c2 + C3) 


2 


10 


10011001 


a(c2 + C3) 


2 



Thus it follows that IC2I = 282, and we constructed code in projective space of length 7 
and minimum distance 3. By adding the whole space and the empty space, we obtain that 
^2(7,3) > 573. 



3.6 Encoding in Grassmannian 

In this section we consider the encoding in Grassmannian. Cover in [20] provided a scheme for 
calculating the index of any sequence in a given subset of binary n-sequences according to its 
position in the lexicographic ordering of this subset. The goal of this section is to find similar 
scheme for subspaces in Grassmannian. The first, intuitive attempt to solve this problem is to 
use Cover's scheme for identifying the vectors of the subspaces, and than define some order on 
the subspaces with the equal identifying vector. But this method is inefficient, since the number 
of subspaces corresponding to each identifying vector is different for each vector, which results 
in inefficient binary representation of indices. 

Therefore, we seek for the one-to-one correspondence between the set of binary vectors of a 
certain length and a set of /c-subspaces, that can be efficient calculated, as has been done for 
other metrics, for example by Knuth in |21j and by Cover in |[20j. 

In the following, we first bound the size of Grassmanian, which is the Gaussian coefficient, 
by some power of 2, in other words, we need to find an integer A'', such that : 



2^ < 



n 
k 



< 2^+\ 



Next, we use the lower bound for finding the correspondence between the set of all binary vectors 
of length N and the subset of the set of all k- subspaces in G2{n, k). Finally, we use the upper 
bound for finding the correspondence between the set of all /c-subspaces in G2{n.,k) and the 
subset of the set of all binary vectors of length -|- 1. 
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Theorem 3.6.1. For 1 < A; < n, 



2k{n—k)+l 



< 



n 
k 



< 



2fc(n-fc)+2 



Proof. The upper bound follows from the Lemma 2.5. In order to prove the lower bound, 

n 
k 



we prove that k{n — k) + 1 < log2{ 

'(2" - l)(2"~i - l)...(2"-'^'+i 



l0g2( 



n 
k 



= log2 



(2-l)(2-2 1)...(2^'-1) 

= log2(2" - 1) + log2(2--i - 1) + ... + log2(2"-'=+i) 

- log2(2 - 1) - log2(22 - 1) - ... - log2(2^ - 1) 

= n- Sn + in-1) - Sn-i + ... + (n - A; + 1) - Sn-k+i 

- O-2 + S2- ...-k + Sk 

= (n + n- 1 + ... + n- A; + 1) - (2 + 3 + ... + A;) 

+ {S2 +S3 + ... + 5k) - {Sn + Sn-1 + ... + Sn-k+l) 

= nk-l-k-{k + l){k - 2) + (^2 +63 + ... + Sk) - {Sn + Sn-1 + ... + Sn-k+l) 

= k{n-k) + l + {S2 + S3 + ... + 5k) - (Sn + Sn-1 + ... + Sn-k+l), 



where Si = log2(2*) - log2(2* - 1). 

We have to prove that < {S2 + S3 + ... + 5k) - {5n + 5n-i + ••. + ^n-fc+i). 
Note that 52 > 5^ > ... > (5„ > 0. We know that 52 ~ 0.4150, ^3 ~ 0.1926, S4 ^ 0.093. 
Therefore, 52 > Sn + 5n-i,S3 > 5n-2, ■■■,5k > 5n-k+i, thus 

< ((52 + (53 + ... + Sk) - {Sn + Sn-1 + ... + (5„_fe+l). 



thus we have proved that k{n — k) + 1 < log2( 



n 
k 



and 



2k{n—k)+l ^ 



n 
k 



< 



2*;(n-fe)+2 



Theorem 3.6.2. (for general q > 2) 



k(n—k) 



< 



Proof. From the Lemma 2.5, for each g. 



k(n—k) 



< 



n 
k 



n 
k 



< 



k{n-k)+l 



< Aq 



k(n—k) 
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Therefore, if q > 4 then 



.k{n—k) 



< 



n 
k 



But if 2 < g < 4 we need another proof. The following proof works for all q > 2: 

n 
k 



We prove that logq{ 



< k{n -k) + l. 



loga( 



n 
k 



log. 



= log,(g^^ - 1) + log,(g"-^ - 1) + ... + log,(g 

- logg(g - 1) - logg(g2 - 1) - ... - log^(g'= - 1) 

= n- 5n + {n-l) - + ... + {n - k + I) - 5n-k^-i 

- 1 + 81-2 + 82- ...-k + 5k 

= (n + n - 1 + ... + n - + 1) - (1 + 2 + 3 + ... + A;) 

+ (81 + 82 + 83 + ... + 8k) - i8n + 8n-l + ... + 8n-k+l) 

= nk - k - k{k - 1) + {81 + 82 + 83 + ... + 8k) - ((5n + 8n-i + ■■■ + 8n-k+i) 

= k{n -k) + {81 + 82 + 83 + ... + 8k) - i8n + 8n-i + ■■■ + 8n-k+l), 



where 8i = logg{q'-) - logq(g* - 1). 
Now we need to prove that 



{81 +82 + 83 + ... + 8k) - {8n + 8n-l + ... + 8n-k+l) < 1, 



but we will prove the stronger claim: 



81 + 82 + 83 + ... + 8k < 1. 



81 + 82 + 83 + ... + 8k = ^log 



i=l 



- 1 



Denote by f{x) = log 



f(x) is monotone decreasing function with lim f{x) = 0. 



Thus we can use the properties of step function and get the following inequality: 

k 



E/(0< / f{x)dx<Y,f{i) 
i=2 •[ i=l 

Denote by g{y) = logy{y^ /y^ — 1), for x > 1 and y > 2. g'{y) < 0, therefore g{y) is monotone 
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decreasing function. Thus, 



oo oo 

I f{x)dx = I 



logo 



7^-1 



dx < J loi 



3^ 



3^-1 



dx 



but we can calculate (using Mathematica) 



oo 

/ 



loga 



3^-1 



dx ^ 0.303421, 



therefore 
k 



< logs 



3-1 



1=2 

CX) 

+ y logs 
1 



3^ - 1 



dx ^ 0.36907 + 0.303421 < 1, 



thus we have proved that log„( 



n 
k 



< k{n — A;) + 1, and 



k{n—k) 



< 



n 
k 



< 



k{n-k)+l 



3.6.1 Encoding in Grassmannian with q=2. 



By Theorem 3.6.1, the longest binary vectors that we can consider are of the length k{n — k) + 1. 
First we consider binary vectors of length k{n — k). 

Let V € GFg*'"' that consists of blocks Vi,V2, ...,Vk of length n — k .Therefore, we can 
write: 

^_ ^ ... ^ 

n—k n—k n—k 

We encode this vector v to the fc-subspace with generator matrix of the following form: 

/ 1 Vi\ 
1 V2 

V 1 Vk J 

Now we want encode all binary vectors of length k{n — A;)+l. First, we recall the following 
Theorem 2.14: 
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Theorem 2.14. |7]. Let 



n 

k 



k(n—k) 



-I 9 



Then the coefficient a/ is the number of partitions of I whose Ferrers diagrams fit a box of 
size khy n — k. 

Corollary 3.6.1.1. The three last coefficients aj^[n-k)i o,k{n-k)-i-, flfc(n-fc)-2 are equal to 
1,1,2, respectively. 

Proof. The only partitions of k{n — k) and k{n — k) — 1 whose Ferrers diagrams fit a box of 
size k X n — k are: 



n-k 



n-k 



Therefore, ak(n-k) = 1 and ak(n-k)-i = 1- 

There are exactly two partitions of ^(n — A;) — 2 whose Ferrers diagrams fits a box of size 
k X n — k: 





n-k 






• 


k 




• 






• 



and 



n-k 



therefore, a^n-k)-2 = 2. 



Now consider the three last summands in the sum of 



n 
k 



in Theorem 2.14: 



IL^ 00^ 

k n — k 

IL^ 01 00^ 

k — 1 n — k — 1 

IL^ 001 00^ 

k-1 n-k -2' 

IL^ Oil 00^ 

k-2 n-k-1 ' 
Thus, we can encode all binary vectors of length k{n — A;) + 1 as follows: 

Let V G GF^'''^'''^^^ , i.e. v = {vi,V2, .■.,v^n-k)+i)- Denote by Vi = {vi, ...,Vn-k), V2 = 

{Vn-k+l, ■■■,V2{n-k)),---yk-l = {'V(k-2){n-k)+li ■■ ■ , V(^k-l){n-k)) , ^k = ('U(fc_l)(n-fc)+l > • • •) «A:(n-fc) ) ■ 

Denote by Vi \ {S} block of size n — k — \S\ that consists of all entries of Vi that are not in S. 



Note, that there are 2^^'^~^'> subspaces with identifying vector 
2k{n-k)-i gubspaces with identifying vector 
2k{n-k)-2 subspaces with identifying vector 
2fe(n-fe)-2 subspaces with identifying vector 
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li V = FiV2"-Vfel, we encode v to the fc-subspace with generator matrix of the following 
form: 

/ 1 \ 
1 V2 

V I Vk J 

If V = ViV2...Vk \ {v^n-k)} lOi W6 encode v to the /c-subspace with generator matrix of 
the following form: 



/ 1 
10 



Vl 
Vn-k+l 








Vi\vi 

V2 \ Vn-k+l 



\ 



1 V(^k_2)(n-k)+l Vk-l\V(k-2){n-k)+l 

V 1 Vk\vk^^_k) J 



li V = ViV2...Vk \ {vk{n-k)j'^k{n-k)-i} 100; cncodc V to the /c-subspace with generator 
matrix of the following form: 



/ 1 
10 



Vl 
Un-k+l 



V2 
Vn-k+2 








Vl \ {vi,V2} 
\ {Vn-k+l,Vn-k+2} 



\ 



1 ■U(fc_2)(n-fe)+l V(fe-2)(n-fc)+2 Vfc_i \ {u(fc_2)(n-fc)+l , ^^(fc-2)(n-fc)+2 } 
yOOOO 1 Vk\{v^n-k),Vk(n-k)-l} J 



li V = ViV2...Vk \ {'yfe(n-fe), 'Ufe(n-fe)-i} 000, we encode v to the /c-subspace with generator 
matrix of the following form: 



/ 1 











Vl 











1 

























1 


^(fc-3)(n-A;)+l 























1 





Vo 

















1 



Vi\vi 

V2 \ Vn-k+l 

Vk-2 \ V(k-3){n-k)+l 
Vk-1 \ V(k-l)(n-k) 



\ 



3.6.2 Spaces encoding in Grassmannian with q=2. 

In the Theorem 3.6.1 we proved that for all 1 < fc < n , 



2k{n—k)+l 



< 



n 
k 



< 



2k{n-k)+2 



Until now we use the left inequality for encoding each vector in GF(2)^^^ to any k- 
subspace of GF(2)". 



38 



Now we want to use all the 



n 
k 



/c-subspaces, or in other words we want to find for each 



J 2 

/c-subspace of GF{2)'^ corresponding vector in GF(2)'^("~'^)+^. 

Every fc-subspace of GF(2)" arises uniquely as the row space of A; x n generator matrix in 
the reduced echelon form over GF{2). For each such matrix there is a binary identifying vector 
of weight k and length n that describes the places of leading entries of the matrix. Now we need 
the following definition: 

Definition 3.6.2.1. For each binary vector X of length n and weight k we define a class of 
reduced echelon form as follows: 

|m G GF(2)^'^" |M in reduced echelon form and X = M^d} . 

In above we proved that if the indices of ones in identifying vector are 1 < < n, 

then the size of the corresponding class of reduced echelon form is 2 2 ^ n-yii+t2+...+iu) 
there are _M^LJi J^fi — ^{-^ _|_ ^2 + ... + ik) arbitrary entries in corresponding matrices. 

Note that by Theorem 2.14 the number of classes of reduced form of size i is the number of 
partitions of i whose Ferrers diagrams fit a box of size k hy n — k. 

The main idea of the following encoding is to find for each matrix M in reduced echelon 
form the corresponding vector V of length k{n — k) + 2 that consists of all arbitrary entries of 
M, and also gives information about the class of reduced form of M. 

We first explicate the idea of the encoding by the following example. 

Example 3.6.2.2. Let n = Q, k = 3, then the correspondence between matrices in reduced 
echelon form and vectors of length k(n — k) + 2 can be summarized in Figure ISTTl 

In general, for subspace with identifying vector of the form 

k n — k 

which in the class of reduced echelon form of size 2^^"'~^\ the corresponding vector is of the 
following form 

and for subspace in class of reduced echelon form of size 2'^^""'^)"' , 1 < ? < k{n — A;), the 
corresponding vector is of the following form: 

XiX2-..Xk(n_k)-iW0Vi-i, 

where Xj, 1 < j < k{n — k) — i, are elements that are in the arbitrary entries of generator 
matrix in echelon form, and Vi-i is the binary vector of the length i — 1 with the following 
properties: 

• Vi-i = 00. ..0, or 

• Vi-i = 00. ..01, or 
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Figure 3.1: Example 3.6.2.2 



matrix in reduced form 



corresponding 
vector 



matrix in reduced form 



corresponding 
vector 



1 Xi X2 

1 X4 X5 xq 
1 Xs Xg 



(xia;2...X7X8a;9lO) 




X2 X3 
1 X5 
1 



(a;ia;2a;3X4 1000000) 



1 xi 2:2 X3 

1 X4 X5 X6 

1 X7 X8 



(X1X2...X6X7X8IOO) 




Xi X2 

1 X3 X4 
1 



(X1X2X3X4IOOOOOI) 



1 Xi X2 X3 
1 X4 X5 X6 

1 X7 



(X1X2...X6X7IOOO) 




Xi X2 
1 X3 
1 X4 



(X1X2X3X41000011) 



1 Xi 






X2 X3 

1 X4 X5 
1 xe X7 



(X1X2...X6X7IOOI) 




Xi X2 
1 X3 

1 



(X1X2X3IOOOOOOO) 



1 Xi X2 X3 

1 X4 X5 xg 
1 



(X1X2...X610000) 




Xi X2 

1 X3 

1 



(X1X2X3IOOOOOOI) 



1 Xl 







X2 X3 

1 X4 X5 

1 X6 



(X1X2...X610001) 




1 Xl 

1 X2 

1 X3 



(X1X2X310000011) 



1 







Xl X2 

1 X3 X4 

1 ,l-5 .T(i 



(X1X2...X6 10011) 




Xl X2 

10 

(J 1 



(xiX2 100000000) 



1 Xl 






X2 X3 

1 X4 X5 
1 



(X1X2X3X4X5 100000) 




1 Xl 

1 X2 

1 



(xiX2 100000001) 



1 Xl 






X2 X3 
1 X4 
1 X5 



(X1X2X3X4X5 100001 




1 Xl 

10 
1 



(xi 1000000000) 



1 





Xl X2 

1 X3 X4 

1 X5 



(X1X2X3X4X5 100011) 




10 
10 
1 



(10000000000) 



• Vi-1 = v\ ^Ivl ^1, where vl ^ = 00.. .0, < ^\ <i-S , and ^ is the vector with no 
two or more successive zeros, < \v2~^\ < i — 3. 

For example, for i = 5 there are only 8 possible such vectors: 

xi...Xfc(„_fc)_5l000000 
xi...Xfe(„_fe)_5l000001 
xi...Xjfc(„_fe)_5 1000011 
xi...Xjfc(„_jfc)_5 1000101 
xi...Xjfc(„_jfc)_5 10001 11 
xi...Xfc(„_fc)_5l001011 
xi...Xfc(„_fc)_5l001101 

Xi...Xfc(„_fc)_5l001111. 
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Note that the number of classes of reduced echelon form of size is upper bounded 

by -P(5) = 7, the unrestricted number of partitions of 5. 

Now the question is if there are enough binary vectors in above form for all subspaces? 

Proposition 3.6.2.3. For all k < n, the number of classes of echelon form of size 
< i < k{n — k), equal the number of classes of echelon form of size 2', < i < k{n — k). 

Proof. By Theorem 2.14, the number of classes of echelon form of size i is the number of 
partitions of i whose Ferrers diagrams fit a box of size k x (n — k). We switch the roles of points 
and empty places in Ferrers diagrams. If there are i points, then there are k(n ~ k) — i empty 
places in Ferrers diagrams. Thus the number of the classes of echelon form of size 2^^'"'^^^^''- is 
equal to number of partitions of i with Ferrers diagrams which fit a box of size k x {n — k). 



Example 3.6.2.4. For k = 3, n = 4, and z = 5 we obtain: 



Ferrers diagrams of partition of 5 in box 
of size 3x4 


Ferrers diagrams of partition of 3*4-5=7 
in box of size 3x4 
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Denote by Vj-i the set of binary vectors Vi-i of length i — 1 with the following properties: 

• Vi-i = 00.. .0, or 

• Vi-i = 00... 01, or 

• Vi-i = v\~^lv2~^l, where vl~^ = 00. ..0, < <i-S , and v^'^ is the vector with no 
two or more consecutive zeros, < \v2~^\ <i — 3. 

Proposition 3.6.2.5. The size of Vi-i is F{i + 1), the (i + l)st Fibonacci number. 

Proof. Denote by /(n) the number of vectors of length n with no two consecutive zeros that 
end by 1. The number of such vectors which start with a " 1" is /(n — 1) and if they start with a 
"0" then they start with "01" and their number is /(n - 2). Thus, f{n) = f{n - 1) + f{n - 2). 
Since /(O) = /(I) = 1, we have that /(n) = F{n+1), the {n+l)st Fibonacci number. Therefore, 

i-2 i-2 i-1 

\Vi-i\ = 1 + ^ f{n) = 1 + ^ F(n + 1) = 1 + ^ F{n) = F{i + 1), 

n=0 n=0 n=l 

where the last equality follows from the following property of Fibonacci numbers: 

F(0) + F{1) + F{2) + ... + F{n) = F{n + 2) - 1. 
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By Theorem 2.14 and Proposition 3.6.2.3, the number of classes of reduced echelon form of 
size 2*''("~'^')^*, < i < k{n — k), is upper bounded by -P(i), the unrestricted number of partitions 
of i. From the Proposition 3.6.2.5, the number of corresponding vectors is F{i + 1). 



The two sequences are: 



i 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 




P{i) 


1 


1 


2 


3 


5 


7 


11 


15 


22 


30 


42 


56 


77 


101 


135 


176 


231 




m 





1 


1 


2 


3 


5 


8 


13 


21 


34 


55 


89 


144 


233 


377 


610 


987 





In order to prove that for each subspace exists corresponding vector of above form, we recall 
the following theorem [8]: 



Theorem 2.11. For all n > 0, the partition function P{n) is less than or equal to the 
(n + l)st Fibonacci number F{n + 1). 

Note that we showed the possible encoding method, but this method is not efficient. To 
device the efficient method we have to define the internal order on the number of partitions of 
a given integer, and also the internal order on the number of corresponding vectors. This is one 
of the issues we will consider in the future work. 

Another problem in this method is that the number of subspaces of the same number m of 
arbitrary entries decreases when m > k{n — k)/2, while the Fibonacci number always increases, 
thus we need only a small amount of corresponding vectors in this case. The possible solution 
for this problem is as follows. 

For given parameters k and n, let x be the minimal number such that 

Q - (ao + ai + . . . + a^_3) < 2^ - (ao2^-2 + ai2^-3 + . . . + 0^.32) , 

where ai, < i < x — 3 are the number of partitions of i whose Ferrers diagrams fit 
a box of size k x n — k, from Theorem 2.14. We use the encoding method above only for 
(ao + oi + . . . + ax-s) subspaces of the largest sizes and the other subspaces will be encoded by 
vectors of length k{n — k) + 2 with x last entries that do not appear in the previous vectors. 

Example 3.6.2.6. Let n = 6 and k = 3. It follows that 

- (1 + 1 + 2) = 2^ - (2^ + 2^ + 2 * 2), 
thus in this case x = 5 and we have the encoding as in Figure 13.21 
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Figure 3.2: Example 3.6.2.6 



matrix in reduced form 



corresponding 
vector 



matrix in reduced form 



corresponding 
vector 



1 ,Ti ,T2 .X3 
1 X4 X5 X6 
1 Xg Xg 



(a;iX2...X7X8X9lO) 




X2 X3 
1 X5 

1 



(X1X2X3X4OOOIIOI) 



1 Xi X2 X3 

1 X4 X5 xq 

1 X7 X8 



(X1X2...X6X7X8100) 




Xi X2 

1 X3 X4 
1 



(X1X2X3X4OOOIIII) 



1 Xi X2 X3 

1 X4 X5 xg 
1 X7 



(X1X2...X6X7IOOO) 




Xi X2 

1 X3 
1 X4 



(X1X2X3X40010000) 



1 Xl 






X2 X3 

1 X4 X5 
1 xe X7 



(X1X2...X6X7IOOI) 




Xl X2 
1 X3 

1 



(X1X2X300010001) 



1 Xl X2 X3 

1 X4 X5 xe 
1 



(X1X2...X600000) 




Xl X2 
1 X3 

1 



(X1X2X300010011) 



1 Xl 






X2 X3 

1 X4 X5 

1 X6 



(X1X2...X6OOOOI) 




1 Xl 

1 X2 
1 X3 



(X1X2X3OOOIOIOI) 



1 






Xl X2 

1 X3 X4 

1 X5 X6 



(X1X2...X6OOOII) 




Xl X2 

10 

1 



(X1X2OOOOIOIII) 



1 Xl 






X2 X3 

1 X4 X5 
1 



(X1X2X3X4X5000101) 




1 Xl 

1 X2 

1 



(xiOOOOOllOll) 



1 Xl 






X2 X3 
1 X4 
1 X5 



(X1X2X3X4X5OOOIII) 




1 Xl 
10 
1 



(xiOOOOOlllOl) 



1 





Xl X2 

1 X3 X4 

1 X5 



(X1X2X3X4X5001011) 




10 

10 
1 



(00000011111) 
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